Marine Envenomation in Okinawa: Overview and Treatment Concept.

Okinawa prefecture is a popular tourist destination due to its beaches and reefs. The reefs host a large variety of animals, including a number of venomous species. Because of the popularity of the reefs and marine activities, people are frequently in close contact with dangerous venomous species and, thus, are exposed to potential envenomation. Commonly encountered venomous animals throughout Okinawa include the invertebrate cone snail, sea urchin, crown-of-thorns starfish, blue-ringed octopus, box jellyfish, and fire coral. The vertebrates include the stonefish, lionfish, sea snake, and moray eel. Treatment for marine envenomation can involve first aid, hot water immersion, antivenom, supportive care, regional anesthesia, and pharmaceutical administration. Information on venomous animals, their toxins, and treatment should be well understood by prehospital care providers and physicians practicing in the prefecture.

New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae).

BACKGROUND: Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates. RESULTS: Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development. CONCLUSIONS: We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

Sexual dimorphism in skull size and shape of Laticauda colubrina (Serpentes: Elapidae).

Background: Sexual dimorphism in size and shape is widespread among squamate reptiles. Sex differences in snake skull size and shape are often accompanied by intersexual feeding niche separation. However, allometric trajectories underlying these differences remain largely unstudied in several lineages. The sea krait Laticauda colubrina (Serpentes: Elapidae) exhibits very clear sexual dimorphism in body size, with previous studies having reported females to be larger and to have a relatively longer and wider head. The two sexes also differ in feeding habits: males tend to prey in shallow water on muraenid eels, whereas females prey in deeper water on congerid eels. Methods: I investigated sexual dimorphism in skull shape and size as well as the pattern of skull growth, to determine whether males and females follow the same ontogenetic trajectories. I studied skull characteristics and body length in 61 male and female sea kraits. Results: The sexes differ in skull shape. Males and females follow distinct allometric trajectories. Structures associated with feeding performance are female-biased, whereas rostral and orbital regions are male-biased. The two sexes differ in allometric trajectories of feeding-related structures (female biased) that correspond to dietary divergence between the sexes. Conclusions: Sea kraits exhibit clear sexual dimorphism in the skull form that may be explained by intersexual differences in the feeding habits as well as reproductive roles. The overall skull growth pattern resembles the typical pattern observed in other tetrapods.

Functional Duplication of the Short-Wavelength-Sensitive Opsin in Sea Snakes: Evidence for Reexpanded Color Sensitivity Following Ancestral Regression.

Color vision is mediated by ancient and spectrally distinct cone opsins. Yet, while there have been multiple losses of opsin genes during the evolution of tetrapods, evidence for opsin gains via functional duplication is extremely scarce. Previous studies have shown that some secondarily marine elapid snakes have acquired expanded "UV-blue" sensitivity via changes at key spectral tuning amino acid sites of the Short-Wavelength Opsin 1 (SWS1) gene. Here, we use elapid reference genomes to show that the molecular origin of this adaptation involved repeated, proximal duplications of the SWS1 gene in the fully marine Hydrophis cyanocinctus. This species possesses four intact SWS1 genes; two of these genes have the ancestral UV sensitivity, and two have a derived sensitivity to the longer wavelengths that dominate marine habitats. We suggest that this remarkable expansion of the opsin repertoire of sea snakes functionally compensates for the ancestral losses of two middle-wavelength opsins in the earliest (dim-light adapted) snakes. This provides a striking contrast to the evolution of opsins during ecological transitions in mammals. Like snakes, early mammals lost two cone photopigments; however, lineages such as bats and cetaceans underwent further opsin losses during their adaptation to dim-light environments.

Comparative Venom Multiomics Reveal the Molecular Mechanisms Driving Adaptation to Diverse Predator-Prey Ecosystems in Closely Related Sea Snakes.

Predator-prey arms races are ideal models for studying the natural selection and adaptive evolution that drive the formation of biological diversity. For venomous snakes, venom is a key bridge linking snakes with their prey, but whether and how venom evolves under the selection of diet remains unclear. Here, we focused on two closely related sea snakes, Hydrophis cyanocinctus and Hydrophis curtus, which show significant differences in prey preferences. Data-independent acquisition (DIA)-based proteomic analysis revealed different degrees of homogeneity in the venom composition of the two snakes, which was consistent with the differential phylogenetic diversity of their prey. By investigating the sequences and structures of three-finger toxins (3FTx), a predominant toxin family in elapid venom, we identified significant differences between the two sea snakes in the binding activity of 3FTx to receptors from different prey populations, which could explain the trophic specialization of H. cyanocinctus. Furthermore, we performed integrated multiomic profiling of the transcriptomes, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and proteomes of the venom glands; constructed venom-related mRNA-miRNA-lncRNA networks; and identified a series of noncoding RNAs involved in the regulation of toxin gene expression in the two species. These findings are highly informative for elucidating the molecular basis and regulatory mechanisms that account for discrepant venom evolution in response to divergent diets in closely related snakes, providing valuable evidence for the study of coselection and coevolution in predator-prey ecosystems.

Identification and Characterization of a Novel Cathelicidin from Hydrophis cyanocinctus with Antimicrobial and Anti-Inflammatory Activity.

The abuse of antibiotics and lack of new antibacterial drugs has led to the emergence of superbugs that raise fears of untreatable infections. The Cathelicidin family of antimicrobial peptide (AMP) with varying antibacterial activities and safety is considered to be a promising alternative to conventional antibiotics. In this study, we investigated a novel Cathelicidin peptide named Hydrostatin-AMP2 from the sea snake Hydrophis cyanocinctus. The peptide was identified based on gene functional annotation of the H. cyanocinctus genome and bioinformatic prediction. Hydrostatin-AMP2 showed excellent antimicrobial activity against both Gram-positive and Gram-negative bacteria, including standard and clinical Ampicillin-resistant strains. The results of the bacterial killing kinetic assay demonstrated that Hydrostatin-AMP2 had faster antimicrobial action than Ampicillin. Meanwhile, Hydrostatin-AMP2 exhibited significant anti-biofilm activity including inhibition and eradication. It also showed a low propensity to induce resistance as well as low cytotoxicity and hemolytic activity. Notably, Hydrostatin-AMP2 apparently decreased the production of pro-inflammatory cytokines in the LPS-induced RAW264.7 cell model. To sum up, these findings indicate that Hydrostatin-AMP2 is a potential peptide candidate for the development of new-generation antimicrobial drugs fighting against antibiotic-resistant bacterial infections.

The banded colour patterns of sea snakes discourage attack by predatory fishes, enabling Batesian mimicry by harmless species.

The evolution of bright 'warning' colours in nontoxic animals often is attributed to mimicry of toxic species, but empirical tests of that hypothesis must overcome the logistical challenge of quantifying differential rates of predation in nature. Populations of a harmless sea snake species (Emydocephalus annulatus) in New Caledonia exhibit colour polymorphism, with around 20% of individuals banded rather than melanic. Stability in that proportion over 20 years has been attributed to Batesian mimicry of deadly snake species by banded morphs of the harmless taxon. This hypothesis requires that banded colours reduce a snake's vulnerability to predation. We tested that idea by pulling flexible snake-shaped models through the water and recording responses by predatory fish. Black and banded lures attracted similar numbers of following fish, but attacks were directed almost exclusively to black lures. Our methods overcome several ambiguities associated with experimental studies on mimicry in terrestrial snakes and support the hypothesis that banded colour patterns reduce a non-venomous marine snake's vulnerability to predation.

Redescriptions of the type specimens of synonymous nominal taxa of sea snakes (Serpentes: Elapidae: Hydrophis, Laticauda) at the Zoological Survey of India.

We redescribe and illustrate the type specimens of ten taxa of sea snakes of the genera Hydrophis Latreille in Sonnini Latreille, 1801 and Laticauda Laurenti, 1768 in the collections of the Zoological Survey of India. The specimens comprise holotypes and syntypes of ten synonymous nominal taxa that represent seven valid nominal taxa. We here clarify that one specimen ZSI 8278 is a syntype of Hydrophis dayanus Stoliczka, 1872, not holotype as previously stated. In one case, four holotypes of four nominal taxa are synonyms of the same taxonHydrophis cyanocinctus Daudin, 1803. Many of these type specimens are herein first depicted in photographs in a publication.

A cathelicidin antimicrobial peptide from Hydrophis cyanocinctus inhibits Zika virus infection by downregulating expression of a viral entry factor.

Zika virus (ZIKV) is a re-emerging flavivirus that causes conditions such as microcephaly and testis damage. The spread of ZIKV has become a major public health concern. Recent studies indicated that antimicrobial peptides are an ideal source for screening antiviral candidates with broad-spectrum antiviral activities, including against ZIKV. We herein found that Hc-CATH, a cathelicidin antimicrobial peptide identified from the sea snake Hydrophis cyanocinctus in our previous work, conferred protection against ZIKV infection in host cells and showed preventative efficacy and therapeutic efficacy in C57BL/6J mice, Ifnar1-/- mice, and pregnant mice. Intriguingly, we revealed that Hc-CATH decreased the susceptibility of host cells to ZIKV by downregulating expression of AXL, a TAM (TYRO3, AXL and MERTK) family kinase receptor that mediates ZIKV infection, and subsequently reversed the negative regulation of AXL on host's type I interferon response. Furthermore, we showed that the cyclo-oxygenase-2/prostaglandin E2/adenylyl cyclase/protein kinase A pathway was involved in Hc-CATH-mediated AXL downregulation, and Hc-CATH in addition directly inactivated ZIKV particles by disrupting viral membrane. Finally, while we found Hc-CATH did not act on the late stage of ZIKV infection, structure-function relationship studies revealed that α-helix and phenylalanine residues are key structural requirements for its protective efficacy against initial ZIKV infection. In summary, we demonstrate that Hc-CATH provides prophylactic and therapeutic efficacy against ZIKV infection via downregulation of AXL, as well as inactivating the virion. Our findings reveal a novel mechanism of cathelicidin against viral infection and highlight the potential of Hc-CATH to prevent and treat ZIKV infection.

Biological and medical aspects related to the yellow-bellied sea snake Hydrophis platurus (Linnaeus, 1766): A view from Colombia.

In Colombia, 317 species of snakes have been recognized, of which 51 (17%) have medical importance due to the toxicity of their venom. A total of 95% of envenomations are caused by snakes of the family Viperidae and 5% of the family Elapidae. The latter form of envenomation is mainly caused by snakes of the genus Micrurus. The only sea snake described is the yellow-bellied snake (Hydrophis platurus), present in the Pacific Ocean. Although Colombia has approximately 1300 km on the Pacific coast and a significant presence of H platurus, envenomation is rare. As a result of the care of a patient with this type of envenomation and of the donation of a H platurus specimen to our laboratory, we decided to conduct this review on the most relevant biological, epidemiological and clinical aspects of this enigmatic and interesting species.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae.

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication-Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

Preparation and Evaluation of a Horse Antiserum against the Venom of Sea Snake Hydrophis curtus from Hainan, China.

Sea snake venom is extremely toxic, and it can induce severe respiratory failure and cause high mortality. The most effective first aid treatment for sea snake bites is to inject antivenom as soon as possible. However, in China, there are only four types of terrestrial snake antivenoms, none of which are effective in the treatment of sea snake bites. In order to develop an antivenom for the dominant species of sea snakes in Chinese seas, Hydrophis curtus venom (HcuV) was chosen as the antigen to immunize horses. From immune plasma, a high-titer Hydrophis curtus antivenom (HcuAV) was prepared. In vitro assessment showed that HcuAV had a cross-neutralizing capacity against HcuV and Hydrophis cyanocinctus venom (HcyV). In vivo assessment indicated that HcuAV injection could significantly improve the survival rates of the HcuV and HcyV envenomated mice (0% to 100% and 87.5%, respectively) when it was injected at a sufficient amount within the shortest possible time. In addition, HcuAV could also effectively alleviate multiple organ injuries caused by HcuV. These results provide experimental support for the future clinical application of HcuAV.

Combined proteomic strategies for in-depth venomic analysis of the beaked sea snake (Hydrophis schistosus) from Songkhla Lake, Thailand.

This study focuses on comprehensive characterization of the venom proteome of the beaked sea snake (Hydrophis schistosus) from Songkhla Lake, Thailand. H. schistosus can be considered as the deadliest sea snake commonly found in the Pacific and Indian oceans. Their envenomation causes muscular paralysis and rhabdomyolysis. To develop effective treatment for this snakebite, it is necessary to understand the detailed venom composition. In this study, multiple mass spectrometry-based approaches were employed. Bottom-up proteomics revealed that tryptic digestion in-solution provided a higher number of toxin proteins identified and a larger sequence coverage, compared to in-gel digestion. In addition, a venom gland transcriptome-derived database was constructed and used as a reference, which 43 known and novel toxin proteins were identified using this database and the UniProtKB. Three-finger toxin and phospholipase A2 were shown to be top two most abundant protein families. Minor compositions included other toxin families and a number of non-toxin proteins. Moreover, a hybrid de novo sequencing was performed to enhance identification of the small proteins/peptides. Using non-digested samples, there were 46 predicted toxin peptides. The finding from this study could lead to a better understanding in pathological effects of the snakebite and the future development of effective antivenoms. SIGNIFICANCE: This study provides a better understanding of the venom proteome composition of the beaked sea snake (H. schistosus) found in the Gulf of Thailand, using a combination of different sample preparation techniques, Serpentes protein database searching, transcriptome-derived protein database searching, and a hybrid de novo peptide sequencing strategy. It revealed 13 toxin protein families and novel proteins in the beaked sea snake venom including new species of phospholipase A2s (PLA2s) and three-finger toxins (3FTxs). It could serve as a basis for the development of snakebite treatments and for the discovery of novel pharmaceutical drugs from the toxin peptides.

Frequency-dependent Batesian mimicry maintains colour polymorphism in a sea snake population.

Evolutionary theory suggests that polymorphic traits can be maintained within a single population only under specific conditions, such as negative frequency-dependent selection or heterozygote advantage. Non-venomous turtle-headed sea snakes (Emydocephalus annulatus) living in shallow bays near Noumea in New Caledonia exhibit three colour morphs: black, black-and-white banded, and an intermediate (grey-banded) morph that darkens with age. We recorded morph frequencies during 18 consecutive years of surveys, and found that the numbers of recruits (neonates plus immigrants) belonging to each morph increased in years when that morph was unusually rare in the population, and decreased when that morph was unusually common. Thus, morph frequencies are maintained by negative frequency-dependent selection. We interpret the situation as Batesian mimicry of highly venomous sea snakes (Aipysurus, Hydrophis, Laticauda) that occur in the same bays, and range in colour from black-and-white banded to grey-banded. Consistent with the idea that mimicry may protect snakes from attack by large fish and sea eagles, behavioural studies have shown that smaller fish species in these bays flee from banded snakes but attack black individuals. As predicted by theory, mimetic (banded) morphs are less common than the cryptically-coloured melanic morph.

Divergence in life-history traits among three adjoining populations of the sea snake Emydocephalus annulatus (Hydrophiinae, Elapidae).

Life-history traits such as rates of growth, survival and reproduction can vary though time within a single population, or through space among populations, due to abiotically-driven changes in resource availability. In terrestrial reptiles, parameters such as temperature and rainfall generate variation in life-histories-but other parameters likely are more important in marine systems. We studied three populations of sea snakes (Emydocephalus annulatus) in adjacent bays in the IndoPacific archipelago of New Caledonia. The extreme philopatry of individual snakes allows us to unambiguously allocate each animal to one of the three populations. Although water temperatures and rainfall do not differ over this small scale, one site experiences more intense winds, restricting opportunities for foraging. Our 18-year mark-recapture dataset (> 1,200 snakes, > 2,400 captures) reveals significant divergence among populations in life-history traits. Survival rates and population densities were similar among sites, but snakes at the most wind-exposed site (Anse Vata) exhibited lower body condition, slower growth, less frequent production of litters, and smaller litters. Weather-driven variation in feeding rates thus may affect life-history traits of marine snakes as well as their terrestrial counterparts, but driven by different parameters (e.g., wind exposure rather than variation in temperatures or rainfall).

Oil spill causes mass mortality of sea snakes in the Gulf of Oman.

Oil spills in the marine environment inflict significant impacts on a wide diversity of marine fauna. Despite the abundance of literature describing these impacts on numerous species, no studies describe the impacts on sea snakes. In this study we report, for the first time, details of an oil spill which caused mass mortality of sea snakes. In this study, 39 sea snake mortalities from the Gulf of Oman, in particular, the coast of Kalba, Sharjah, UAE, were examined. The investigated sea snakes belong to four different species (Hydrophis platurus, H. lapemoides, H. spiralis and H. ornatus). The majority (84.6%) of sea snakes were observed to have oil covering 75-100% of their bodies. The majority (91.4%) of sea snakes were also observed with oil covering their snouts and eyes. A large proportion (25.8, 41.4 and 34.5%) of sea snakes were observed with oil in their mouth, esophagus and stomach.

Sexual dimorphism in size and shape of the head in the sea snake Emydocephalus annulatus (Hydrophiinae, Elapidae).

In snakes, divergence in head size between the sexes has been interpreted as an adaptation to intersexual niche divergence. By overcoming gape-limitation, a larger head enables snakes of one sex to ingest larger prey items. Under this hypothesis, we do not expect a species that consumes only tiny prey items to exhibit sex differences in relative head size, or to show empirical links between relative head size and fitness-relevant traits such as growth and fecundity. Our field studies on the sea snake Emydocephalus annulatus falsify these predictions. Although these snakes feed exclusively on fish eggs, the heads of female snakes are longer and wider than those of males at the same body length. Individuals with wider heads grew more rapidly, reproduced more often, and produced larger litters. Thus, head shape can affect fitness and can diverge between the sexes even without gape-limitation. Head size and shape may facilitate other aspects of feeding (such as the ability to scrape eggs off coral) and locomotion (hydrodynamics); and a smaller head may advantage the sex that is more mobile, and that obtains its prey in narrow crevices rather than in more exposed situations (i.e., males).

Population dynamics of the sea snake Emydocephalus annulatus (Elapidae, Hydrophiinae).

For sea snakes as for many types of animals, long-term studies on population biology are rare and hence, we do not understand the degree to which annual variation in population sizes is driven by density-dependent regulation versus by stochastic abiotic factors. We monitored three populations of turtle-headed sea snakes (Emydocephalus annulatus) in New Caledonia over an 18-year period. Annual recruitment (% change in numbers) showed negative density-dependence: that is, recruitment increased when population densities were low, and decreased when densities were high. Windy weather during winter increased survival of neonates, perhaps by shielding them from predation; but those same weather conditions reduced body condition and the reproductive output of adult snakes. The role for density-dependence in annual dynamics of these populations is consistent with the slow, K-selected life-history attributes of the species; and the influence of weather conditions on reproductive output suggests that females adjust their allocation to reproduction based on food availability during vitellogenesis.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon.

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 µg/g) and long-chain neurotoxin (i.p. 0.14 µg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Venom-gland transcriptomic, venomic, and antivenomic profiles of the spine-bellied sea snake (Hydrophis curtus) from the South China Sea.

BACKGROUND: A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes. RESULTS: A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A2 (PLA2), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA2 (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA2s. CONCLUSIONS: Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.